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We analyze the running gauge coupling at finite temperature for QCD, using the functional 
renormalization group. The running of the coupling is calculated for all scales and temperatures. 
At finite temperature, the coupling is governed by a fixed point of the 3-dimensional theory for 
scales smaller than the corresponding temperature. The running coupling can drive the quark 
sector to criticality, resulting in chiral symmetry breaking. Our results provide for a quantitative 
determination of the phase boundary in the plane of temperature and number of massless flavors. 
Using the experimental value of the coupling at the r mass scale as the only input parameter, we 
obtain, e.g., for Nt = 3 massless flavors a critical temperature of T cr ss 148 MeV in good agreement 
with lattice simulations. 
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I. INTRODUCTION 


Strongly interacting matter is believed to have funda¬ 
mentally different properties at high temperature than at 
low or zero temperature lj. Whereas the latter can be 
described in terms of ordinary hadronic states, a hadronic 
picture at increasing temperature is eventually bound to 
fail; instead, a description in terms of quarks and gluons 
is expected to arise naturally owing to asymptotic free¬ 
dom. In the transition region between these asymptotic 
descriptions, effective degrees of freedom, such as order 
parameters for the chiral or deconfining phase transition, 
may characterize the physical properties in simple terms, 
i.e., with a simple effective action 0. 

If a simple description at or above the phase transi¬ 
tion does not exist and the system is strongly interacting 
in all conceivable sets of variables [||, a formulation in 
terms of microscopic degrees of freedom has the greatest 
potential to bridge wide ranges in parameter space from 
first principles. 

In this Letter, we report a nonperturbative study of 
finite-temperature QCD parameterized in terms of mi¬ 
croscopic degrees of freedom: gluons and quarks. We 
use the functional renormalization group (RG) [1, [B, f| 
and concentrate on two problems which are accessible in 
microscopic language: first, we compute the running of 
the gauge coupling driven by quantum as well as ther¬ 
mal fluctuations, generalizing previous zero-temperature 
studies 0. Second, we investigate the induced quark dy¬ 
namics including its back-reactions on gluodynamics, in 
order to monitor the status of chiral symmetry at finite 
temperature. With this strategy, the critical tempera¬ 
ture of chiral symmetry restoration can be computed. 
This Letter is particularly devoted to a presentation of 
the central physical mechanisms which our study reveals 
for the dynamics near the phase transition; all technical 
details can be found in a separate publication @]. 

The functional RG yields a flow equation for the effec¬ 


tive average action T*, Q, 

d t T k = ^STrd t R k (T^+R k )~\ t = ln|, (1) 

where T^ interpolates between the bare action T k=A = S 

( 2 ) 

and the full quantum effective action T = Tz—o; T k 
denotes the second functional derivative with respect to 
the fluctuating field. The regulator function Rk specifies 
the Wilsonian momentum-shell integration, such that the 
flow of Tfc is dominated by fluctuations with momenta 
p 2 ~ k 2 . 

An approximate solution to the flow equation can reli¬ 
ably describe also nonperturbative physics if the relevant 
degrees of freedom in the form of RG relevant operators 
are kept in the ansatz for the effective action. As the cru¬ 
cial ingredient, the choice of this ansatz has to be guided 
by all available physical information. 

II. TRUNCATED RG FLOW FOR THERMAL 
GLUODYNAMICS 

In this work, we truncate the space of possible ac¬ 
tion functionals to a tractable set of operators which is 
motivated from various sources and principles. For the 
principle of gauge invariance, we use the background- 
field formalism as developed in [9j, i.e., we work in the 
Landau-de Witt background-field gauge and follow the 
strategy of 0,0 for an approximate resolution of the 
gauge constraints a Decomposing the gauge field into 
a background-field part and a fluctuation field, this strat¬ 
egy focusses on the flow in the background-field sector of 
the action and neglects an independent running of the 
fluctuation-field sector. In general, the solution of the 
strongly-coupled gauge sector represents the most deli¬ 
cate part of this study, owing to a lack of sufficient a 
priori control of nonperturbative truncation schemes. A 
first and highly nontrivial check of a solution is already 
given by the stability of its RG flow, since oversimplify¬ 
ing truncations which miss the right degrees of freedom 
generically exhibit IR instabilities of Landau-pole type. 
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The IR stability of our solution arises from an impor¬ 
tant conceptual ingredient: we optimize our truncated 
flow with an adjustment of the regulator to the spec¬ 
tral flow of r( 2 ) [H|, instead of a naive canonical 

momentum-shell regularization. For this, we integrate 
over shells of eigenvalues of I^ 2 ), by inserting r( 2 ) into 
the regulator and accounting for the flow of these eigen¬ 
values. More precisely, we use the exponential regula¬ 
tor i of the form R k ( F^) = r( 2 )/[exp(r( 2 )/-Z'fc& 2 ) — 1], 
where Z k denotes the wave function renormalization of 
the corresponding field (gluons or ghost). In a perturba¬ 
tive language, the optimizing spectral adjustment allows 
for a resummation of a larger class of diagrams. 

The main part of our truncation consists of an infi¬ 
nite set of operators given by powers of the Yang-Mills 
Lagrangian, 

i\ym[A] = [ e = \f; v f; v . (2) 

J x ^ 

In the function 144(0) = W\B + \W 2 O 2 + ^W^O 3, ..., the 
coefficients IF; form an infinite set of generalized cou¬ 
plings. This truncation represents a gauge-covariant gra¬ 
dient expansion in the field strength, neglecting higher- 
derivative terms and more complicated color and Lorentz 
structures. Hence, the truncation includes arbitrar¬ 
ily high gluonic correlators projected onto their small- 
momentum limit and onto the particular color and 
Lorentz structure arising from powers of F 2 . Since per¬ 
turbative gluons are certainly not the true degrees of free¬ 
dom in the IR, an inclusion of infinitely many gluonic 
operators appears mandatory in order to have a chance 
to capture the relevant physics in this gluonic language. 
The covariant gradient expansion does not only facilitate 
a systematic classification of the gluonic operators, it is 
also a consistent expansion in the framework of the func¬ 
tional RG. 1 

Furthermore, our truncation includes standard (bare) 
gauge-fixing and ghost terms, neglecting any nontrivial 
running in these sectors. We emphasize that we do not 
expect that this truncation reflects the true behavior in 
these sectors, but we assume that the non-trivial run¬ 
ning in these sectors does not qualitatively modify the 
running of the background-held sector where we read off 
the physics. 

A well-known problem of gauge-covariant gradient ex¬ 
pansions in gluodynamics is the appearance of an IR un¬ 
stable Nielsen-Olesen mode in the spectrum [lS]. At fi¬ 
nite temperature T, this problem is severe, since such 
a mode will be strongly populated by thermal fluctua¬ 
tions, typically spoiling perturbative computations fl4| . 
Our how equation allows us to resolve this problem with 


the aid of the IR regulator. We remove this mode’s un¬ 
physical thermal population by a T-dependent regulator. 
In this way, we obtain a strictly positive spectrum for the 
thermal huctuations. 

In the present truncation, the how equation results in 
a differential equation for the function 144 of the form 

d t W k (6) = F[dgW k , d 2 W k , dtdeWk, d t d 2 e W k ], (3) 


where the extensive functional T depends on derivatives 
of 144, on the coupling g and the temperature T; it is dis¬ 
played in Q. We use the nonrenormalization of the prod¬ 
uct of coupling and background held, gA : for a nonper- 
turbative definition of the running coupling in terms of 
the background wave function renormalization Z k = IFi 


/V = d t g 2 = rig 2 , rj =--^-d t Z k . (4) 

£ k 

The how of Z kl and thus the running of the coupling, 
is successively driven by all generalized couplings IF,;. 
Keeping track of all contributions from the hows of the 
IF,, Eq. © boils down to a recursive relation, 


d t W i = f ij {g,T)d t W h (5) 

with fij(g,T) representing the expansion coefficients of 
the RHS of Eq. ([3]), which obey f,j = 0 for j > i + 1. 
Solving Eq. (0 for d t Z k = d t W\, we obtain a nonper- 
turbative (3 g z function in terms of an inhnite asymptotic 
but resumable power series, 

00 ( 2\m+l 

P 9 2 = H a ™(f) r 2 (47r) 2 jm’ ( 6 ) 

with temperature-dependent coefficients a m 2 For ex¬ 
plicit representations of the a m and further details, we 
refer the reader to [§]. At zero T, the /3 g 2 function 
agrees well with perturbation theory for small coupling, 
reproducing one-loop exactly and two-loop within a few- 
percent error. For larger coupling, the resumed integral 
representation of Eq. © reveals a second zero of the f3 g 2 
function for hnite g 2 , corresponding to an IR attractive 
non-Gaufiian hxed point g 2 > 0, which conhrms the re¬ 
sults of [?}]. The resulting zero-temperature how of the 
coupling is displayed by the (red) solid “T=0” line in 
Fig. [T] Deviations from perturbation theory become sig¬ 
nificant at scales below 1 GeV; the IR hxed-point be¬ 
havior sets in on scales of order 0(100 MeV). As initial 
condition, we use the measured value of the coupling at 
the r mass scale [n|, ct s = 0.322, which evolves to the 
world average of a s at the Z mass scale. We stress that 


1 The gluonic gradient expansion as a local expansion can, of 
course, not be expected to give reliable answers to all questions; 
for instance, bound-state phenomena such as glue balls are en¬ 
coded in the nonlocal pole structure of higher-order vertices. 


2 An iteration of the same procedure would result in (3 functions 

for the other higher-order couplings Wi with i > 1. We neglect 

these, since they do not exert a direct influence on the quark 

sector discussed below. 
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FIG. 1: Running SU(3) Yang-Mills coupling ctYM{k,T) as a 
function of k for T = 0, 100, 300 MeV compared to the one- 
loop running for vanishing temperature. 


no other parameter or scale is used as an input neither 
at T = 0 nor at finite temperature as described below. 

The appearance of an IR fixed point in Yang-Mills the¬ 
ories is a well-investigated phenomenon in the Landau 
gauge E3, where the running coupling can be defined 
with the aid of the ghost-gluon vertex. Whereas the 
universal perturbative running of the coupling is iden¬ 
tical for the different gauges despite the different defini¬ 
tions of the coupling, we moreover observe a qualitative 
agreement between the Landau gauge and the Landau- 
de Witt background-field gauge in the nonperturbative 
IR in the form of an attractive fixed point. This points 
to a deeper connection between the two gauges which 
deserves further study and may be traced back to cer¬ 
tain non-renormalization properties in the two gauges. 
Note that the IR fixed point in the Landau gauge is in 
accordance with the Kugo-Ojima and Gribov-Zwanziger 
confinement scenarios [181 ] . In general, an IR fixed point 
is also compatible with the existence of a mass gap 0, 
since such a gap in the physical spectrum typically in¬ 
duces threshold and decoupling behavior towards the IR. 

At finite temperature T, the UV behavior remains un¬ 
affected for scales k T and agrees well with the one- 
loop perturbative running coupling at zero temperature, 
as expected. 3 In the IR, the running is strongly mod- 


3 Our truncation does not reproduce higher orders in the high- 

temperature small-coupling expansion which proceeds with odd 
powers in g beyond one loop Il9ll . These odd powers are a result 
of a resummation which, in the language of the effective action, 
requires non-local operators, being neglected so far. In any case, 
we do not expect the underlying quasi-particle picture of these 
operators to hold near the chiral phase transition, such that their 
omission should not qualitatively modify our low-temperature 
results. 


ified: The coupling increases towards lower scales until 
it develops a maximum near k ~ T. Below, the cou¬ 
pling decreases according to a power law g 2 ~ k/T, 
see Fig. [T] This behavior has a simple explanation: 
the wavelength of fluctuations with momenta p 2 < T 2 
is larger than the extent of the compactified Euclidean 
time direction. Hence, these modes become effectively 
3-dimensional and their limiting behavior is governed by 
the spatial 3d Yang-Mills theory. This dimensional re¬ 
duction has been discussed for the running coupling in 
a perturbative weak-coupling framework in [20| . Our re¬ 
sults generalize this to arbitrary couplings. As a nontriv¬ 
ial new result, we observe the existence of a non-Gaufiian 
IR fixed point g 2 d * also in the reduced 3-dimensional the¬ 
ory. By virtue of a straightforward matching between the 
4d and 3d coupling, the observed power law for the 4d 
coupling is a direct consequence of the strong-coupling IR 
behavior in the 3d theory, g 2 (k <C T) ~ g 2 dif k/T. We 
find a 3d fixed-point value of a^d,* = — 2.7 

which demonstrates that the system is strongly coupled 
despite the naive decrease of the 4 d coupling; also the 3d 
background anomalous dimension is large, approaching 
r] 3 d ■ 1 near the IR fixed point. This scenario is remi¬ 
niscent to hot 4 d </> 4 theory which approaches a strongly 
coupled IR limit at high temperatures analogous to the 
3d Wilson-Fisher fixed point [2l[. The observation of an 
IR fixed point in the 3d theor y ag ain agrees with recent 
results in the Landau gauge j22[. Note that our flow 
to the 3d theory is driven by thermal as well as quan¬ 
tum fluctuations. This is different from a purely ther¬ 
mal flow [H, [24j | as used in (25|, where the IR limit is 
four-dimensional, being characterized by a decoupling of 
fluctuation modes owing to thermal masses. 

The 3d IR fixed point and the perturbative UV be¬ 
havior already qualitatively determine the momentum 
asymptotics of the running coupling. Phenomenologi¬ 
cally, the behavior of the coupling in the transition re¬ 
gion at mid-momenta is most important, which is quan¬ 
titatively provided by the full 4d finite-temperature flow 
equation. 


III. TRUNCATED RG FLOW FOR CHIRAL 
QUARK DYNAMICS 

Extending our calculations to QCD, we first include 
the quark contributions to all gluonic operators of our 
truncation, as done in Ref. (2|| for QED. This effectively 
corresponds to Heisenberg-Euler-type quark-loop contri¬ 
butions to the flow of the function 144(0). Successively, 
we obtain quark-loop contributions to the coefficients a m 
in Eq. ([6]) and thus to the running coupling, accounting 
for the screening nature of fermionic fluctuations; here, 
we confine ourselves to massless quarks, but current- 
quark masses can straightforwardly be included f||. 

The determination of the critical temperature T cr 
above which chiral symmetry is restored requires a sec¬ 
ond crucial ingredient for our truncation: we study the 
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gluon-induced quark self-interactions of the type 

r^int ^ X a ^§1p 0 ,'lljp1p^'lps^ (7) 

where a, /3,... denote collective indices including color, 
flavor, and Dirac structures. The resulting flow equa¬ 
tions for the A’s are a straightforward finite-temperature 
generalization of those derived and analyzed in [H], HH 
and will not be displayed here for brevity. The boundary 
condition X a p-ys —> 0 for k —> A —> oo guarantees that the 
A’s at k < A are solely generated by quark-gluon dynam¬ 
ics from first principles (e.g., by 1PI “box” diagrams with 
2-gluon exchange). We emphasize that this is an impor¬ 
tant difference to, e.g., the Nambu-Jona-Lasinio (NJL) 
model, where the A’s are independent input parameters. 

We consider all linearly-independent four-quark inter¬ 
actions permitted by gauge and chiral symmetry. A pri¬ 
ori, these include color and flavor singlets and octets in 
the (S— P), (V—A) and (V+A) channels. UA(l)-violating 
interactions are neglected, since they may become rele¬ 
vant only inside the xSB regime or for small Nf. We drop 
any nontrivial momentum dependencies of the A’s and 
study these couplings in the point-like limit A(|pi| -C k). 
This is a severe approximation, since it inhibits a study 
of QCD properties in the chirally broken regime; for in¬ 
stance, mesons manifest themselves as momentum singu¬ 
larities in the A’s. Nevertheless, the point-like truncation 
can be a reasonable approximation in the chirally sym¬ 
metric regime, as has recently been quantitatively con¬ 
firmed for the zero-temperature chiral phase transition 
in many-flavor QCD [27j • Our truncation is based on the 
assumption that quark dynamics both near the finite-T 
phase boundary as well as near the many-flavor phase 
boundary [29? | is driven by similar mechanisms. Our re¬ 
strictions on the four-quark interactions result in a total 
number of four linearly-independent A couplings; all oth¬ 
ers channels are related to this minimal basis by means 
of Fierz transformations. 

Introducing the dimensionless couplings A = k 2 A, the 
(3 functions for the A couplings are of the form 

d t X = 2A - AAA - bXg 2 - c 5 4 , (8) 

where the coefficients A, &, c are temperature dependent, 
A being a matrix and b a vector in the s pac e of A cou¬ 
plings (for explicit representations, see [§. I30|). 

Within this truncation, a simple picture for the chiral 
dynamics arises: at weak gauge coupling, the RG flow 
generates quark self-interactions of order A ~ g A via the 
last term in Eq. © with a negligible back-reaction on 
the gluonic RG flow. If the gauge coupling in the IR 
remains smaller than a critical value g < g CI , the A self¬ 
interactions remain bounded, approaching fixed points A* 
in the IR. Technically, the ~ g 4 term is balanced by the 
first term ~ 2A at these fixed points. The fixed points are 
the counter-parts of the GauBian fixed point A® aufi = 0 
in NJL-like models (at g 2 = 0), here being modified by 
the gauge dynamics. At these fixed points, the fermionic 
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FIG. 2: Running QCD coupling a s (k,T) for Nf = 3 massless 
quark flavors and N c = 3 colors and the critical value of the 
running coupling a CI (k, T) as a function of k for T = 130 MeV 
(upper panel) and T = 220 MeV (lower panel). The existence 
of the (a s ,Qcr) intersection point (marked by a circle) in the 
former indicates that the \'SB quark dynamics can become 
critical. 


subsystem remains in the chirally invariant phase which 
is indeed realized at high temperatures T > T cr . 

If the gauge coupling increases beyond the critical cou¬ 
pling g > g CI , the IR fixed points A* are destabilized and 
the quark self-interactions become critical min. Then 
the gauge-fluctuation-induced A’s have become strong 
enough to contribute as relevant operators to the RG 
flow, with the term ~ AAA dominating Eq. ([8|). In this 
case, the A’s increase rapidly, approaching a divergence 
at a finite scale k = k x sb- In fact, this seeming Landau- 
pole behavior indicates ,\SB and the formation of chiral 
condensates: the A’s are proportional to the inverse mass 
parameter of a Ginzburg-Landau effective potential for 
the order parameter in a (partially) bosonized formula¬ 
tion, A ~ 1 /to 2 [HI, t22|- Thus, the scale at which the 
self-interactions formally diverge is a good measure for 
the scale k x sb where the effective potential for the chiral 
order parameter becomes flat and is about to develop a 
nonzero vacuum expectation value. 

Whether or not chiral symmetry is preserved by the 
ground state therefore depends on the running QCD cou¬ 
pling g relative to the critical coupling g cr which is re¬ 
quired to trigger xSB. For instance, at zero tempera¬ 
ture, the SU(3) critical coupling for the quark system is 
a cr = <?c r /(4 tt) ~ 0.8 in our RG scheme [32j, being only 
weakly dependent on the number of flavors [27j . Since 
the IR fixed point for the gauge coupling is much larger 
a* > a cr for not too many massless flavors, the QCD vac¬ 
uum is characterized by ySB. At finite temperature, the 
running of the gauge coupling is considerably modified in 
the IR. Moreover, the critical coupling is T dependent, 
g cr = g C r{T/k). This can be understood from the fact 
that all quark modes acquire thermal masses and, thus, 
stronger interactions are required to excite critical quark 
dynamics. This thermal decoupling is visible in the coef- 
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ficients A , b , and c in Eq. ([8|>. all of which vanish in the 
limit T/k —> co. 

In Fig. [2j we show the running coupling a s and its 
critical value a CI for T = 130 MeV and T = 220 MeV 
as a function of the regulator scale k. The intersection 
point k CI between both marks the scale where the quark 
dynamics become critical. Below the scale k cr , the sys¬ 
tem runs quickly into the ySB regime. We estimate the 
critical temperature T cr as the lowest temperature for 
which no intersection point between a s and a CI occurs. 4 
Compared to Q, we have further resolved the finite-T 
Lorentz structure of the four-fermion couplings (301], re¬ 
sulting in a slightly improved estimate for T cr : we find 
T cr « 172+ 4 ° MeV for N f = 2 and T cr « 148+^ MeV for 
Nf = 3 massless quark flavors in good agreement with 
lattice simulations [H]. The errors arise from the ex¬ 
perimental uncertainties on a s fl6 ]. Dimensionless ratios 
of observables are less contaminated by this uncertainty 
of a s . For instance, the relative difference for T CI for 

lVf=2 and 3 flavors is ~ atI- 3 — = 0.150 .. .0.165 

(T„ f - 2 +T„ f - 3 )/2 

in reasonable agreement with the lattice value 5 of ~ 
0.121 ± 0.069. 

Furthermore, we compute the critical temperature for 
the case of many massless quark flavors Nf, see Fig. [3] 
We observe an almost linear decrease of the critical 
temperature for increasing Nf with a slope of AT cr = 
T(Nf) — T(Nf + 1) ss 24 MeV for small Nf. In addition, 
we find a critical number of quark flavors, Nf 1 = 12, 
above which no chiral phase transition occurs. This re¬ 
sult for Nf r agrees with other studies based on the 2-loop 
/3 function [29|; however, the precise value of iV f cr is ex¬ 
ceptionally sensitive to the 3-loop coefficient which can 
bring Nf* down to V f cr ~ 10 [27|. Since we do not 

consider our truncation to be sufficiently accurate for a 
precise estimate of this coefficient, our study does not 
contribute to a reduction of the current error on iV f cr . 
Instead, we would like to emphasize that the flatten¬ 
ing shape of the phase boundary near iV f cr is a generic 
prediction of the IR fixed-point scenario: here, the sym¬ 
metry status of the system is governed by the fixed-point 
regime where dimensionful scales such as Aqqd lose their 
importance (§j. In any case, since V f cr is smaller than 
V f af ' = V-jV c = 16.5, our study provides further evi¬ 
dence for the existence of a regime where QCD is chiral 
symmetric but is still asymptotically free. 


4 Strictly speaking, this is a sufficient but not a necessary criterion 
for chiral-symmetry restoration. In this sense, our estimate for 
T cr is an upper bound for the true T CI . Small corrections to this 
estimate could arise, if the quark dynamics becomes uncritical 
again by a strong decrease of the gauge coupling towards the IR. 

5 The large uncertainty on the lattice value arises from the fact 
that the statistical errors on the Nf = 2 and Nf = 3 results for 
T cr are uncorrelated. 



FIG. 3: Chiral-phase-transition temperature T CI versus the 
number of massless quark flavors Nf. In the dashed-line re¬ 
gion, we expect UA(l)-violating operators to become quan¬ 
titatively important. The flattening at Nf > 10 is a conse¬ 
quence of the IR fixed-point structure Q. 

IV. CONCLUSION 

In summary, we have determined the ySB phase 
boundary in QCD in the plane of temperature and flavor 
number. Our quantitative results are in accord with lat¬ 
tice simulations for Nf = 2, 3. For larger Nf, we observe 
a linear decrease of T cr , leveling off near V f cr owing to the 
IR fixed-point structure of QCD. Our results are based 
on a consistent operator expansion of the QCD effective 
action that can systematically be generalized to higher 
orders. 

The qualitative validity and the quantitative conver¬ 
gence of this expansion are naturally difficult to analyze 
in this strongly-coupled gauge system, particularly for 
the gluonic sector. The fact that our truncation results in 
a stable RG flow at strong interactions is already a highly 
non-trivial check that any ansatz which misses the true 
degrees of freedom generically fails. A more quantitative 
evaluation of the validity of our expansion will require 
the inclusion of higher-order operators in the covariant 
gradient expansion as well as higher-order ghost terms. 
An inclusion of operators that distinguish between elec¬ 
tric and magnetic sectors at finite T, e.g., (u^F ^) 2 with 
the heat-bath four-velocity u^, should facilitate to dis¬ 
tinguish between differing coupling strengths in the two 
sectors, as done in (2f| using an ansatz inspired by hard 
thermal loop computations. 

We observe an improved control over the truncation in 
the quark sector at least for the chirally symmetric phase, 
which suffices to trace out the phase boundary. Quanti¬ 
tatively, this has been confirmed by a stability analysis 
of universal quantities such as Nf 1 under a variation of 
the regulator in (27J which gives strong support to the 
point-like truncation of the quark self-interactions. Qual¬ 
itatively, the reliability of the quark truncation can also 
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be understood by the fact that the feed-back of higher- 
order operators, such as ~ (V’t/’) 4 , is generally suppressed 
by the one-loop structure of the flow equation. 

Future extensions should include mesonic operators 
which can be treated by RG rebosonization techniques 
& This would not only provide access to the broken 
phase and mesonic properties, but also permit a study 
of the order of the phase transition. For further phe¬ 


nomenology, the present quantitative results that rely on 
only one physical input parameter can serve as a promis¬ 
ing starting point. 
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